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The low-lying isovector dipole strengths in neutron rich nuclei 26 Ne and 28 Ne are 

Hi ■ 

' investigated in the quasiparticle relativistic random phase approximation. Nuclear 

ground state properties are calculated in an extended relativistic mean-field theory 

>. 

, plus BCS method where the contribution of the resonant continuum to pairing corre- 



lations is properly treated. Numerical calculations are tested in the case of isovector 
dipole and isoscalar quadrupole modes in the neutron rich nucleus 22 O. It is found 
that in present calculation low-lying isovector dipole strengths at E x < 10 MeV in 
nuclei 26 Ne and 28 Ne exhaust about 4.9% and 5.8% of the Thomas-Reiche-Kuhn 
dipole sum rule, respectively. The centroid energy of the low-lying dipole excitation 



o 
in 
o 

22 



X 



is located at 8.3 MeV in 26 Ne and 7.9 MeV in 28 Ne. 
PACS numbers: 21.60.Jz, 24.30. Cz, 24.30.Gd 

I. INTRODUCTION 

Nuclear giant resonances have been known since 50 years ago for the dipole mode 
and more than 30 years for the other modes. But the research field was limited to the 
excitations of ouolei aion g the ^ability Hnefl fl Q. Recent,, the radlo act,ve ion 
beam physics has become one of the frontiers in nuclear physics. It offers the possibility 
to broaden the study of the giant resonance to weakly bound nuclei. Nuclei close to 
the drip line present some unique properties: a small separation energy of the valence 
nucleon, the smearing density distribution and a strong coupling between the bound 
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state and the particle continuum. Those exotic properties attract more attentions both 
experimentally and theoretically. Low-lying electric dipole modes may appear in these 
weakly bound nuclei, which are so-called Pygmy Dipole Resonances. Although carrying 
only a small fraction of the full dipole strength these states are of particular interest 
because they are expected to reflect the motion of the neutron skin against the core 
formed with an equal number of protons and neutrons. Recent experiments have shown 
that the increase of the dipole strength at low energies in neutron-rich nuclei could affect 
the corresponding radiative neutron capture cross section considerably^], which has a 
significance in astrophysics. Over the last decade, much experimental and theoretical 
efforts have been dedicated to investigate properties of the low-lying dipole mode in light 
neutron-rich nuclei, in particular, to answer whether or not these dipole excitations can 



be attributed to the collectivit 
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Recently, Beaumel et al.[17j have measured the inelastic scattering of 26 Ne + 208 Pb 
using a 60 MeV/u 26 Ne secondary beam at RIKEN. This reaction is dominated by 
Coulomb excitations and selective for El transitions. The experimental data are now 
under analysis jl^]. As a subsequent work, they will continue the experiment using a more 
neutron rich projectile 28 Ne. Therefore the theoretical investigation of low-lying dipole 
modes in 26 Ne and 28 Ne has a practical significance. 26 Ne and 28 Ne are neutron rich nuclei, 
whose Fermi surfaces are close to the particle continuum. Therefore the description of 
those nuclei has to explicitly include the coupling between bound states and the particle 



continuum. The contribution of the particle continuum to the nuclear 



properties at 
22, Q- On 



low-energies can mainly be attributed to a few resonant states jiol I201 121 
the other hand, it is well known that pairing correlations play an important role in 
describing properties of open shell nuclei. In order to depict the collective excitations of 
those nuclei pairing correlations have to be taken into account. Recently, a number of 
theoretical works has devoted to study the properties of low-lying dipole modes in the 



24, 



22]- 



framework of the quasiparticle random phase approximation (QRPA) 
Paar and his co-workers [2^] have studied the evolution of the low-lying dipole strength 
in Sn-isotopes in the quasiparticle relativistic random phase approximation (QRRPA) in 
the configuration space formalism. 

In this paper, we aim at the investigation of the properties of low-lying dipole modes 
in neutron rich nuclei 26 Ne and 28 Ne in the QRRPA which is formulated in the response 



function method. The QRRPA is an extension of the fully consistent RRPA 



tie response 
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by taking into account the effect of pairing correlations. A consistent treatment of RRPA 
within the RMF approximation requires the configurations including not only the pairs 
formed from the occupied Fermi states and unoccupied states but also the pairs formed 
from the Dirac states and occupied Fermi states. It has been emphasized in Refs. 3, 0] 
that the inclusion of configurations built from the positive energy states in the Fermi sea 
and negative energy states in the Dirac sea is essential to give correctly a quantitative 
description in the excitation energies of isoscalar giant multipole resonances as well as to 
ensure the current conservation and decouple the spurious states. In present calculations, 
we pay more attentions to the energy weighted moment ni\ and the centroid energy of the 
low-lying dipole strength as well as the contribution of states around the Fermi surface 



to the low-lying dipole strength. Although some theoretica 



investigations show that 



26 Ne and 28 Ne are deformed and strongly anharmonic [32|, |33j, |3J] , a spherical symmetry 
is assumed in the present investigation. In order to show the applicability of the method 
with a spherical assumption we also study the quadrupole excitations in these nuclei and 
compare the calculated position and transition strength of the lowest 2 + states with the 
experimental data. 

In this work, the ground state properties of neutron rich nuclei 26 Ne and 28 Ne are calcu- 
lated in the extended relativistic mean field and Bardeen-Cooper-Schrieffer (RMF+BCS) 



approxiamtion 



221 ] . where the resonant continuum is properly treated. The empirical 



pairing gaps deduced from odd-even mass differences are adopted in the BCS calculation 
in this work. All calculations are performed with the parameter set NL3 3fil |. which gives 
a good description of not onl y the g round state properties but also the collective 



giant resonance 



29. 



30, 



31 



37, 



38, 



The paper is arranged as follows. In Sec. II we present the formalism of the QRRPA 
in the extended RMF+BCS ground state in the response function approach. A test of 
the numerical calculation of the QRRPA in the neutron rich nucleus 22 O is performed 
and compared with the results in Refj3], which is given in Sec. III. In Sec. IV the ground 
state properties of 26 Ne and 28 Ne are studied in the extended RMF+BCS approach. 
Then the QRRPA in the response function formalism is applied to study the properties 
of isovector giant dipole resonances in nuclei 26 Ne and 28 Ne. Finally we give a brief 
summary in Sec.V. 
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II. THE QUASIPARTICLE RELATIVISTIC RANDOM PHASE 

APPROXIMATION 

There are usually two methods to obtain the RPA strength in the study of nuclear 

collective excitations. One is working in a particle-hole configuration space and solving 

I — I 

the RPA equation by a matrix diagonalization method 30], and another is based on 



the linear response theory |40j]. In the response function formalism one solves a Bethe- 
Salpeter equation by inversion. In both methods the starting point is a self-consistent 
solution of the nuclear ground state. In this paper we shall work in the response function 
formalism and study nuclear dipole excitations in neutron rich nuclei. 

In the RRPA calculation we first solve the Dirac equation and equations of meson 
fields self-consistently in the coordinate space. The continuum is discretized by ex- 
panding nucleon spinors on a complete set basis, such as eigenfunctions in a spherical 
harmonic oscillator potential. Those single particle states are used to build the RRPA 
configurations: a set of particle-hole pairs(p/i) and pairs formed from the negative energy 
state in the Dirac sea and the hole state in the Fermi sea(a/i). 

The response function of a quantum system to an external field is given by the imag- 
inary part of the polarization operator: 

R(Q, Q- k, k'; E) = -ImU R (Q, Q- k, k'; E) , (1) 

7T 

where Q is an external field operator. The RRPA polarization operator is obtained by 
solving the Bethe-Salpeter equation: 

n(Q,Q;k,k',E)= U (Q,Q;k,k',E)-J29i [ d 3 k 1 d 3 k 2 n. (Q, r*; k, ki, E) 

A(ki, k 2 , E)U(Ti, Q; k 2 , k', E) , (2) 

In the RRPA, the residual particle-hole interactions are obtained from the same La- 
grangian as in the description of the nuclear ground state. They are generated by ex- 
changing various mesons: the isoscalar scalar meson a, the isoscalar vector meson uj 
and the isovector vector meson p. Therefore, in Eq.(2) the sum i runs over a, uo and p 
mesons, and and Di are the corresponding coupling constants and meson propagators. 
They are T* = 1 for a meson and T l = 7^,7^3 for uj and p mesons, respectively. The 
meson propagators in the non-linear model are non-local in the momentum space, and 
therefore have to be calculated numerically. The detailed expressions of non-linear meson 



propagators A(ki, k 2 , E) can be found in Ref.J^. n is the unperturbed polarization 
operator, which in a spectral representation has the following retarded form: 



Il*(P,Q;k,k';E) 




h.a 



(MPlHpWQlHh) (<p p \\P L H h )(4> h \\QLH P ) 

E-{e p - e h ) +irj E + (e p — e h ) + it] 
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E 



[Sa ~ Eh) + iV 



E 



[Sa - Eh) + ir] 



(3) 



The unperturbed polarization operator includes not only the positive energy particle- 
hole pairs but also pairs formed from the Dirac sea states and Fermi sea states. In 
RefsflQ 

the authors show that the inclusion of configurations built from the positive 
energy states in the Fermi sea and negative energy states in the Dirac sea is essential 
to give correctly a quantitative description in the excitation energies of isoscalar giant 
multipole resonances as well as to ensure the current conservation and decouple the 
spurious states. 

The pairing correlation and coupling to the continuum are important for exotic nuclei. 
A proper treatment of the resonant continuum to pairing correlations has been rece ntly 
investigated in the Hartree-Fock (HF) Bogoliubov or the HF+BCS approximation 1 1 9. 
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2jj and the extended RMF+BCS approximation [22, 23]. It shows that the simple 



BCS approximation in the resonant continuum with a proper boundary condition works 
well in the description of ground state properties even for neutron rich nuclei. We shall 
treat the pairing correlation in the BCS approximation in this work and the resonant 
continuum is calculated by imposing an asymptotic scattering boundary condition. 

When pairing correlations are taken into account, the elementary excitation is a two- 
quasiparticle excitation, rather than a particle-hole excitation. The unperturbed polar- 
ization operator in the QRRPA in the response function formalism can be constructed 
in a similar way: 



n*(P,Q;k,k';E) 



' o,/9 



M\PLUp)(<l>p\\QLUa) (tpWPLMaMJQLHp) 



a,P 



E — (E a + Ep) +i V E + (E a + E p ) + ir] 

(UPd^^QLUh) {^\\P L Ua)m\QL\\<t>p)' 



E-{E a + \-£j 



IT] 



E + (E a + A - eg) + ir] 



(,4) 
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with 

A aP = (u a v p + {-l) L v a up) 2 {l + 5 a p)~ l , (5) 

where v 2 is the occupation probability and u 2 a = 1 — v 2 . E a = a/ (s a — A) 2 + A 2 is 
the quasiparticle energy, where A and A are the Fermi energy and pairing correlation 
gap, respectively. In the BCS approximation, the 4> a is the eigenfunction of the single 
particle Hamiltonian with an eigenvalue e a . In Eq.(4), terms in the first square bracket 
represent those excitations with one quasiparticle in fully or partial occupied states and 
one quasiparticle in partial occupied or unoccupied states. Terms in the second square 
bracket describe all excitations between positive energy fully or partial occupied states 
and negative energy states in the Dirac sea. For unoccupied positive energy states 
outside the pairing active space, their energies are Ep = e@ — A, occupation probabilities 
Vp = and it| = 1. For fully occupied positive energy states, the quasiparticle energy 
and the occupation probability are E a = A — e a and v 2 = 1 in Eq.(4). States in the 
Dirac sea are not involved in pairing correlations. Therefore those quantities and vL 
are set to be and 1, respectively. Once the unperturbed polarization operator in the 
quasiparticle scheme is built, the QRRPA response function can be obtained by solving 
the Bethe-Salpeter equation (2) as usually done in the RRPA. 



III. NUMERICAL CALCULATION AND TEST OF THE QRRPA 

In this section, we first check the validity of the present QRRPA calculations. We 
apply the QRRPA to calculate the response function of the isovector giant dipole res- 
onance(IVGDR) and the isoscalar giant quadrupole resonance(ISGQR) in the neutron 
rich nucleus 22 . Similar calculations for the nucleus 22 were recently performed by 
Paar et aljssj in the framework of the Relativistic Hartree-Bogoliubov(RHB) + QRRPA 
in the configuration space formalism. 

The ground state properties of the nucleus 22 are calculated in the extended 



RMF+BCS approach|22j with the parameter set NL3. The neutron pairing gap is ob- 
tained from the experimental binding energies of neighboring nuclei, A n = 1.532 MeV. 
In the QRRPA calculation particle-hole residual interactions are taken from the same ef- 
fective interaction NL3, which is used in the description of the ground state of 22 O. Fully 
occupied states and states in the pairing active space are calculated self-consistently in 
the extended RMF+BCS approach in the coordinate space. The BCS active space is 
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FIG. 1: IVGDR strengths in 22 O. The solid curve represents the result calculated in the 
QRRPA approach. The result performed in the RRPA approach is shown by a dashed curve. 
All results are calculated with the effective Lagrangian parameter set NL3. 

taken as all states in the sd shell as well as U7/2 state, which is a resonant state in 22 O. 
A scattering boundary condition is imposed in the resonant continuum. Unoccupied 
states outside of the pairing active space are obtained by solving the Dirac equation in 
the expansion on a set of the harmonic oscillator basis. The response functions of the 
nuclear system to the external operator are calculated at the limit of zero momentum 
transfer. It is also necessary to include the space-like parts of vector mesons in the QR- 



4l|. The consistent 



RPA calculations, although they do not play role in the ground state 
treatment guarantees the conservation of the vector current. 

In Fig.l, we show the response function of the IVGDR mode in 22 O calculated in the 
RRPA and QRRPA approaches. The isovector dipole operator used in the calculations 

is3 : 



N Z Z N 

Q = e— - ^2 r i Y iM(ri) - e— riY 1M (fi) 



(6) 

i=i ~ i=i 

which excites an L = 1 type electric (spin-non-flip) AT = 1 and AS = giant resonance 

with J w = 1~. The spurious state for exotic nuclei may appear at the energy around 1 

MeV in the IVGDR strength in numerical calculations due to the mixture of the isoscalar 



mode[43[. In our present calculations, the spurious state is removed by slightly adjusting 
the coupling constant of the a meson in the residual interaction by less than 1%, that 
does not affect the general results. 

In general, the IVGDR strengths in light stable nuclei are expected to be fragmented 
substantially. This also occurs in the response function of the IVGDR in neutron rich 
nuclei. More fragmented distributions around the GDR region in 22 are observed in 
Fig.l. In addition to the characteristic peak of the IVGDR at the energy around 20 MeV, 
the low-lying dipole strength appears at the excitation energy below 10 MeV. It can be 
seen that the inclusion of pairing correlations enhances the low-lying dipole strength 
and has a slight effect on the strength at the normal dipole resonance. This illustrates 
the importance of including pairing correlations in the study of the low-lying isovector 
dipole strength in neutron rich nuclei. The effect of pairing correlations on the isovector 
dipole strength in 22 observed in our calculation is consistent with that obtained in the 
RHB+QRRPA in the configuration formalism (Fig.2 of Ref.Q). 
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FIG. 2: ISGQR strengths in 22 0. Notations are the same as in Fig.l. 



The response functions of the ISGQR mode in 22 calculated in the RRPA and 
QRRPA approaches are shown in Fig.2. In present calculation the isoscalar quadrupole 
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TABLE I: The energy weighted moment mi at E x < 60 MeV for the electric isovector dipole 
and isoscalar quadrupole excitations in 22 0. DC represents the result from double commutators 
in the non-relativistic approach. 





DC 




mi(E x <60 MeV) 






RRPA 


QRRPA 


IVGDR(e 2 .fm 2 MeV) 


75.9 


82.08 


81.71 


ISGQR(e 2 /m 4 MeV) 


2018 


1971 


2159 



operator is taken from Ref. 



42] 



(7) 



The inclusion of pairing correlations shifts the low-lying quadrupole strength to higher 
energy region and enhances the low-lying quadrupole strength, while this only slightly 
affect the strength at the normal giant resonance region. A similar result on the isoscalar 
quadrupole strength in 22 O has been observed in Fig. 3 of Ref. [3] in the framework of 
the RHB+QRRPA. 

In Table I we show the energy weighted moment mi at E x <60 MeV for electric 
isovector dipole and isoscalar quadrupole excitations in 22 0. DC represents the result 
from double commutators in the non-relativistic approach. In the isovector dipole mode 
the value corresponds to the Thomas-Reiche-Kuhn (TRK) dipole sum rule. It is shown 
that the mi obtained by the integration of the RPA strength till to 60 MeV is slightly 
larger than that obtained by the double commutator (DC) in the dipole mode. While in 
the quadrupole mode both RRPA and QRRPA results are close to the DC value. 



IV. ISOVECTOR DIPOLE EXCITATION IN NEUTRON RICH NUCLEI 

26 NE AND 28 NE 

A. Ground state properties of nuclei 26 Ne and 28 Ne 

Ground state properties of nuclei 26 Ne and 28 Ne are studied in the extended 
RMF+BCS with the parameter set NL3, where a spherical symmetry is assumed. The 
continuum is calculated by imposing a scattering boundary condition and the width of 
the resonant state is not considered in this work. Constant pairing gaps are adopted in 
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TABLE II 
properties: 
rms radii, 
energy 44 ] . 



Neutron and proton pairing gaps in 26 Ne and 28 Ne, and calculated ground state 
neutron and proton Fermi energies, binding energies as well as neutron and proton 
Values in the parenthesis are the corresponding experimental data of the binding 



20 



Ne 



28 



Nc 



A n (MeV) 
A p (MeV) 
Xn(MeV) 
X p (MeV) 
E B {MeV) 
r n (fm) 



r p (fm) 



1.436 
2.025 
-5.325 
-14.168 
201.8(201.6) 
3.179 
2.784 



1.400 
2.101 
-4.290 
-16.247 
210.5(206.9) 
3.348 
2.833 



the calculation of pairing correlations, which are obtained from the experimental binding 
energies of neighboring nuclei by the formula: 

A p = - (B(Z - 2, N) - AB(Z - 1, N) + 6B{Z, N) - AB{Z + 1,N) + B{Z + 2, N)) , 



A n = - (B(Z, N - 2) - 4B(Z, N - 1) + 6B(Z, N) - 4B(Z, N + 1) + B(Z, N + 2)) . 
8 

(9) 

In our calculations, the neutron pairing active space in nuclei 26 Ne and 28 Ne includes 
states up to the iV = 28 major shell and 2p 3 / 2 state, which are ld 5 / 2 , 2si/ 2 , ld 3 /2, 2p 3 / 2 , 
and 1/7/2- The BCS active space for proton is taken as all states in the sd shell. In 
Table II we list the neutron and proton pairing gaps in nuclei 26 Ne and 28 Ne derived 
from Eqs.(8,9) and the calculated ground state properties, including neutron and proton 
Fermi energies, total binding energies as well as the neutron and proton rms radii. The 
values in the parenthesis are corresponding experimental binding energies taken from 



Ref. 



44 1 . Neutron single particle energies and BCS occupation probabilities for those 



states near the neutron Fermi energy are shown in Table III, where levels (2p 3 / 2 and 
I/7/2) with positive energies are the single particle resonant states. 



11 



TABLE III: Neutron single-particle energies e a and occupation probabilities -u 2 of levels near 
the neutron Fermi energy in nuclei 26 Ne and 28 Ne. 





26 Ne 

e a (MeV) 


< 

ex. 


28 Ne 

e a (MeV) 


vl 

Ot 


14/2 


-10.548 


0.982 


-10.836 


0.989 


2«l/2 


-6.549 


0.824 


-7.054 


0.946 


14/2 


-3.408 


0.099 


-4.299 


0.503 


2j>3/2 






0.786 


0.018 


1/7/2 


2.946 


0.007 


2.223 


0.011 




r(fm) 

FIG. 3: Neutron and proton density distributions in nuclei 26 Ne and 28 Ne. All results are 
calculated in the RMF+BCS. 

Thus the nucleon density with pairing correlations can be written as: 

p(r)=J2^^v 2 Mr)Mr), (10) 

a 

where the summation runs over all states weighted by the factor v 2 a . In Fig. 3 we show 
the calculated nucleon densities in 26 Ne and 28 Ne. It is shown in Table II that proton rms 
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TABLE IV: The calculated energies and B(E2) values for the lowest 2 + states in 26 Ne and 28 Ne. 
The values in the parenthesis are the corresponding experimental data taken from Ref . . 





E 2 +(MeV) 


B(E2){e 2 fm A ) 


26 Ne 


1.46(1. 99±0.012) 


223(228±41) 


28 Ne 


1.29(1.32±0.020) 


156(269±136) 



radii are much smaller than those of neutron. The neutron densities have far extending 
tails seen in Fig. 3, which clearly shows neutron skins formed in those nuclei. 

B. Low-lying isovector dipole modes in 26 Ne and 28 Ne 

In this work a spherical assumption is adopted in the RMF and RRPA calculation, 
although 26 Ne and 28 Ne are deformed and anharmonic. In order to see how far the 
present model can be effective one first studies the quadrupole excitations in these nuclei 
and compare with the experimental data. The calculated energies and B(E2) values for 
the lowest 2 + states in 26 Ne and 28 Ne are listed in Table IV. The values in the parenthesis 



are the corresponding experimental data taken from Ref. 45[. Although the B(E2) value 
in 28 Ne is very close to the lower limit of experimental data, the present calculations 
reasonably reproduce the lowest 2 + states and its B(E2) values. 

We now apply the QRRPA approach to investigate the isovector dipole response in 
26 Ne and 28 Ne. We focus our attention on properties of the isovector low-lying dipole 
strength. In Fig.4 we present the Hartree and perturbed strengths for the isovector dipole 
mode in nuclei 26 Ne (upper panel) and 28 Ne (lower panel) calculated in the RRPA and 
QRRPA approaches. Short-dashed curves represent the RRPA strengths, the QRRPA 
response functions are denoted by solid ones. In Fig.4, it can be seen that the perturbed 
strengths at the energy above 10 MeV for those neutron rich nuclei are very fragmented. 
Compared to the Hartree strengths, the RPA strengths are shifted to higher energy region 
due to the repulsive particle-hole residual interaction generated mainly by exchanging p 
meson. 

In addition to the characteristic peak of the IVGDR around energy of 20 MeV, low- 
lying dipole strengths appear at the excitation energy below 10 MeV. It is shown in 
Fig.4 that the effect of pairing correlations on the isovector dipole strength in 26 Ne 
and 28 Ne shifts the low-lying dipole strength to higher energy region and decreases the 
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low-lying dipole strength, especially in 26 Ne. The situation slightly differs from that in 
22 0, where the low-lying strength is increased. The decrease of the low-lying strength is 
mainly due to the fact that two-quasiparticle excitations in the dipole mode are weakened 
by a factor of v 2 a when the pairing correlation is switched on. In addition, the two- 
quasiparticle excitation en ergy is larger than the corresponding particle-hole excitation 
energy in the dipole mode 27[ . On the contrary, when the pairing correlation is taken 
into account the configuration space becomes larger, which allows for the particle-particle 
and hole-hole transitions. This enlarged configuration space may increase the low-lying 
strength. As analyzed below, the state 2si/ 2 produces a large contribution to the low- 
lying strength of the dipole mode in 26 ' 28 Ne and 22 0. It is known that the 2sx/2 state is 
located below and above the Fermi surface in 26 ' 28 Ne and 22 0, respectively. Therefore 
particle-particle excitations, especially due to the particle state 2si/ 2 , largely enhance 
the low-lying strength in 22 0, which is not true in the case of 26 Ne and 28 Ne. 

In comparison with the Hartree strength it is found that the RPA strength of the 
isovector dipole mode at the low energy region in 26 Ne and 28 Ne shown in Fig.4 remains 
at its position and is contributed mainly from a few particle-hole configurations, which 
shows a single particle like property. The low-lying strength is slightly attracted back to 
the lower energy, which is due to the correlations of the isoscalar operator in the isovector 
mode 4^ . Differing from the normal IVGDR response, the low-lying resonance can be 
interpreted as the excitation of the excess neutrons out of phase with the core formed 
with an equal number of protons and neutronsQ]. Analyzing the Hartree strength of 
the isovector dipole mode in 26 Ne at the energy below 10 MeV calculated with pairing 
correlations, one finds a pronounced peak around 8.5 MeV, which is formed mainly from 
the neutron configurations of u(2s~l2p 3/2 ) (8.482 MeV) and u(2s~l2p l/2 ) (9.232 MeV), 
where the value in the parenthesis is its Hartree energy. Since 16 neutrons in 26 Ne fill 
neutron orbits up to 2s 1/2 and form a sub-closed shell, the occupation probabilities at 
lc?3/2 and I/7/2 states in 26 Ne are relatively small, see Table III. Therefore the contribu- 
tion from neutron states 1^3/2 and I/7/2 to the low-lying Hartree strength is insignificant. 
In contrast the Hartree strength at the low energy region in 28 Ne, in addition to the peaks 
formed from the neutron configurations of u(2sy 2 2p 3 / 2 ) (8.523 MeV) and i/(2s^ / 1 2 2p 1 / 2 ) 
(9.080 MeV), a few more peaks appear, which are formed from the neutron excitation 
between the bound level ld^/ 2 and levels in the continuum. It is found that the Hartree 
strengths at the low-lying dipole in 26 Ne and 28 Ne are mainly due to neutron excitations 
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FIG. 4: Isovector dipole strength functions in neutron rich nuclei 26 Ne (upper panel) and 
28 Ne (lower panel). The QRRPA responses with the pairing (solid curves) are compared with 
the RRPA calculation without the pairing (short-dashed curves). The thick and thin curves 
represent perturbed and Hartree strengths, respectively. 
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FIG. 5: The contribution of proton to the isovector dipole strength in 26 Ne (left panel) and 
28 Ne (right panel) in the QRRPA approach. The solid curves are the full QRRPA strengths. 
The short-dashed curves represent the strengths from proton. 

TABLE V: The non-energy weighted moment mo and energy weighted moment m\ of isovector 
dipole strengths in 26 Ne and 28 Ne in the QRRPA calculations. We separate the energy region 
into the low energy (0 MeV < E x < 10 MeV) and the high energy (10 MeV < E x < 30 
MeV). The values in the last two column are obtained from the classical TRK dipole sum rule 
and TRK cluster sum rule(e 2 /m 2 MeV). The units are e 2 fm 2 and e 2 /m 2 MeV for mo and mi, 
respectively. 





0MeV< E x 


<10MeV 


10MeV< E x 


<30MeV 








m 


m\ 


mo 


m± 


Strk 


Scius 


26 Ne 


0.542 


4.525 


5.032 


103.9 


91.7 


17.2 


28 Ne 


0.705 


5.606 


5.648 


111.2 


95.8 


21.3 



near the Fermi surface. 

Although the Hartree low-lying strength is mainly formed from the neutron excita- 
tions, the RPA strengths are fully correlated and contributed from both neutron and 
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TABLE VI: Centroid energies of the isovector dipole response functions in 6 Ne and °Ne. The 



centroid energies are calculated within ~ 


1U MeV ana 1U 


~ 60 MeV, respectively. All energy 


values are in unit of MeV. 








26 Ne 


28 Ne 


E(0 ~ 10) 


8.34 


7.94 


£(10 ~ 60) 


22.32 


21.13 



proton. To illustrate the contribution of the proton to the full strength we set a very 
small value of the neutron effective charge in the dipole operator instead of eZ/A. Re- 
sults are plotted in Fig.5. The short-dashed curves represent the strengths from proton, 
which are compared with the full QRRPA strengths denoted by solid curves. It clearly 
shows that the proton also plays an important role in the perturbed strength even at 
low-lying dipole states. 



In Ref. 



14j . the authors have studied the evolution of collectivity in the isovector 



dipole response in the low- lying region for neutron-rich isotopes of O, Ca, Ni, Zr, and Sn 
using RRPA method. They conclude that in light neutron-rich nuclei, such as neutron- 
rich isotopes of O and Ca, the onset of dipole strength in the low-lying region is due 
to single-particle excitations of the loosely bound neutrons. By analyzing the structure 
of RRPA strengths in low-lying region in nuclei 26 Ne and 28 Ne, we find that there are 
only several configurations contributing to the low-lying RRPA strengths and the RRPA 
strengths remain their positions compared to the Hartree strengths in the low-energy 
region, which mean that the RRPA strengths in low-lying region in nuclei 26 Ne and 28 Ne 
are dominated by single-particle transitions. 

In order to give a more clear description of those isovector low-lying dipole states 
obtained in the QRRPA approach, we calculate various moments of isovector dipole 
strengths at a given energy interval: 

m k = [ 2 R L (E)E k dE , (11) 



where E\ and E 2 are the lower and upper energies of the integral, respectively. The RPA 
equation is solved till E = 60 MeV in the present calculations. The non-energy weighted 
moment m and the energy weighted moment m\ of isovector dipole strengths in 26 Ne 
and 28 Ne are calculated at two energy intervals: the low energy region (0 MeV < E x < 
10 MeV) and the high energy region (10 MeV < E x < 30 MeV), which are listed in Table 



17 



IV. The values obtained from the Thomas- Reiche-Kuhn (TRK) dipole sum rule are listed 
in the sixth column. In our present QRRPA calculations the low-lying isovector dipole 
strengths in 26 Ne and 28 Ne exhaust about 4.93% and 5.85% of the TRK dipole sum rule, 
respectively. This is consistent with the recent experimental observations in 18 O, 20 O 
and 22 O, where the low-lying isovector dipole strengths exhaust about 5% of the TRK 
dipole sum rulejl^j]. In general the percentage of the low- lying isovector dipole strength 
becomes large as the increase of the neutron excess jl^. The energy weighted moment 
m 1 at E x < 10 MeV in 28 Ne is about 1.0% larger than that in 26 Ne. 

On th e other hand, th e cluS te r SUm rdefl^HQ - U8Ually US ed to .nde.W the 
properties of low-lying dipole strength in exotic nucleus excitations. Here we choose 20 Ne 
as the core. In Table V we also list the values obtained form TRK cluster dipole sum rule, 
say 17.2 e 2 /m 2 MeV and 21.3 e 2 /m 2 MeV for 26 Ne and 28 Ne, respectively. It shows that 
the low-lying dipole excitations in these two neutron rich nuclei exhaust about 26.3% of 
the TRK cluster dipole sum rule. Similar results are obtained in Ref . [z| for neutron rich 
Oxygen isotopes. 

The centroid energy of the response function is defined as: 

E = m 1 /m . (12) 

We separate the energy interval into two regions: MeV< E x <10 MeV and 10 MeV< 
E x <60 MeV. The centroid energies in these two energy regions are listed in Table V. In 
the QRRPA calculations the centroid energies of low-lying isovector dipole strengths in 
26 Ne and 28 Ne are 8.34 MeV and 7.94 MeV, respectively. Whereas the centroid energies 
of the normal IVGDR strengths are located at 22.3 MeV for 26 Ne and 21.1 MeV for 28 Ne. 



V. SUMMARY 



In this paper we have studied the properties of low-lying isovector dipole resonances 
in neutron rich nuclei 26 Ne and 28 Ne in the framework of the QRRPA with the effec- 
tive Lagrangian parameter set NL3. The ground state properties are calculated in the 
extended RMF+BCS approach, where the resonant continuum is properly treated. Con- 
stant pairing gaps extracted from the experimental binding energies of neighboring nuclei 
are adopted in the BCS calculation. In the QRRPA calculation the negative energy states 
in the Dirac sea are included due to the completeness. It is shown that the inclusion 
of pairing correlations has a relatively strong effect on the low-lying isovector dipole 
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strength in neutron rich nuclei. In the QRRPA calculation the low-lying isovector dipole 
strengths in 26 Ne and 28 Ne exhaust about 5% and 26.3% of the TRK dipole sum rule 
and TRK dipole cluster sum rule, respectively. The centroid energies of the low-lying 
dipole excitation in nuclei 26 Ne and 28 Ne are located at the energy around 8.0 MeV. 
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